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Abstract 

We study the origin of non-analyticity in a s of a short-distance QCD observable to 
demonstrate that the infrared renormalons, the same-sign factorial growth of the per- 
turbative expansion, is a universal phenomenon that originates entirely from the small 
coupling domain. In particular, both the position and the nature of the singularity of 
the Borel transform of the perturbative series prove to be independent of whether the 
running coupling a(k 2 ) becomes singular at some finite scale ("Landau pole"), or stays 
finite down to k 2 = 0. We argue that getting hold of the infrared renormalons per se can 
help next to nothing in quantifying non-perturbative effects. 
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1 Introduction 



Last years showed revival of interest to the problem of asymptotic behaviour of the perturbative 
(PT) expansion in QCD |L|. Considering the Borel transform of PT series has been suggested 
as technical means of improving convergence, and demonstrated that the PT expansion by itself 
prompts about nontrivial dynamics at low momentum scales via infrared (IR) renormalons. 

In reality, one does not expect @ QCD observables to have proper analyticity in a s in the vicinity 
of a s = 0, which is necessary to justify the Borel resummation as a mathematically well defined 
operation. The original arguments rely on the nontrivial analytic properties of QCD observables 
with respect to the external momentum scale parameter Q 2 . These properties are driven by causality 
and unitarity constraints, supplemented with the pattern of the hadron spectrum. On the other 
hand, the IR renormalons themselves have a simpler origin related merely to the nature of the 
perturbative expansion. In the present paper we address these features of the PT expansion to 
illustrate that the IR renormalons are not a reflection of the actual non-perturbative (NP) dynamics. 
Rather they are an artefact of an attempt to describe physics occurring at quite different scales by 
means of one and the same expansion parameter. The salient conclusions we draw from our analysis 
are: 

• IR renormalons persist even if the coupling stays finite at arbitrarily small momentum scales 
when no apparent uncertainty associated with the Landau singularity shows up. 

• The position and the nature of the IR renormalon singularity in the Borel image of a generic 
hard QCD observable are determined by the first two coefficients of the /3-function. 

• The modification of the Borel integration prescription one can design to represent the true 
answer, depends heavily on the details of strong interaction dynamics as well as on the particu- 
lar observable. Therefore, such a possibility seems to carry no practical value as a perturbative 
algorithm. 

• The Q 2 -dependence of a short distance observable, inferred from the Borel resummation, 
typically yields an incorrect image of the actual magnitude of condensate effects, when considered 
at intermediate Q 2 . 

We also briefly address the OPE-motivated procedure to illustrate that it is free from the IR 
renormalon problems. 

In this paper we outline the main framework and the results of the analysis supplying mini- 
mal illustrations when necessary. A more complete discussion with better elaborated qualitative 
and quantitative considerations and deeper physical arguments about the relevance of the analysis 
undertaken, will be presented in a forthcoming publication ||. 

2 Toy model 

To study the origin of the IR renormalons and their relation to NP dynamics we concentrate on 
a simplified model peeled off inessential details. To this end, we address the problem of the PT 
expansion of the integral 




as a model for a short-distance dominated observable. Integrals of this type naturally emerge, for 
example, when one calculates the one-gluon correction to a collinear safe QCD observable using the 
running coupling in the integrand. Correctness of this substitution can be unambiguously proven 
in an Abelian theory; in the QCD context such a doing is often motivated by the "naive non- 
Abelization" or by the "extended BLM prescription" for fixing the relevant hardness scale @ in 
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Feynman integrals. 

It is easy to see that the observable (2A) satisfies the inhomogeneous renormalization group 
(RG) equation 

+ P J («) = a, (2.2) 



where 

da(k 2 ) 



(3 (a(k 2 )) . (2.3) 



d log k 2 

Hereafter, to simplify notation, we absorb the factor fto/Air into the coupling and denote by a 
(without an argument) its value at the external hard scale of the problem, 

Arc 

In this notation the PT expansion of the QCD /^-function is 

(3(a) = —a 2 — 7a 3 + higher terms , 7 = ^ > . (2.4) 

Po 

In this Section we concentrate on the properties of / as a function of a (its perturbative expansion, 
Borel representation, analyticity in a). The analytic properties of 1(a) stem from the dependence 
of a(k 2 ) on a at k 2 < Q 2 , which dependence is determined by the RG trajectories 

«(* 2 ) da' Q 2 _ 

In what follows we shall refer to t as "time" . 

It is clear that the integral fl2.1|) is sensible only if a(k 2 ) does not develop a Landau singularity 
at finite positive k 2 (at finite time), which implies that (3(a) /a has a zero on the positive real axis. 
(The first such zero will be generically denoted by a, < a < +00.) Arguments in favour of an 
"infrared finite" coupling as the only reasonable expansion parameter for QCD observables, at least 
in the present context, will be given in a more detailed publication |3[]. 

Our aim is to compare the "exact" expression Q2.1p with the results one obtains using the Borel 
resummation tricks. 

2.1 Analyticity in a: perturbative and physical "phases" 



It is straightforward to solve the RG equation (|2.2|) . First, one finds the solution of the homogeneous 
equation 

dAj>) + = Q = 

da 

Function f(a) is nothing but the RG- invariant expression for a pure power of the momentum scale. 
For small a, for example, one invokes the expansion ( |2.4| ) to obtain, up to an overall factor, 

f(a) = exp {-1} cT* (1 + O (a)) = . (2.7) 
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In terms of f(a) one can write the solution of the inhomogeneous equation in the form 

r a da' n/ 1 

1(a) = f p (a) -^-^rry (2.8) 
J aoo p{a!) f p {a') 



Here a^ is the asymptotic value of the running coupling at t = oo along the trajectory ( |2.5|) that 
starts from a given a at t — 0. We shall call such point(s) attractive. In the physical phase, that is 
for real positive (small) initial a values, a^ = a : 

rPH , . f a da' a' [ f a da" ) ,„ . 

' (a) = Lw) ew r p L^n\- (2 - 9) 



From explicit expressions (|2.8| ), ( |2.7|) it is clear that an analytic continuation of 1(a), starting from 
some small positive a, can fail - for small a - only due to non-analyticity in a of the lower limit 
of the integral, that is a OC3 (a). Let us demonstrate that it does fail; the faster, the smaller initial 
value of a is taken. 

To this end let us look at the RG trajectories a(t) for complex initial values of a. We start with 
the case when \a\ is small but its phase is finite. In this case the value of \a(t)\ stays uniformly 
small along the whole trajectory, 0<t<oo, so that it suffices to approximate the full /3-function by 
its one-loop expression to find 

- -t)[l + 0(a(t)\na(t))] « ~-t. (2.10) 



a(t) \a J a 

In this approximation trajectories lie on small circles (either in the upper or in the lower half-plane) 
that touch the real axis at a = 0. The radius of a circle is related to the phase of the initial a by 



Ima 



-i 



i-i^L < l. (2.11) 
2 Ima V ; 



It is important to realize that the existence of the family of such trajectories having the common 
limiting point a^ = d = — is a universal, purely perturbative feature of the asymptotically free 
theory with (3 oc —a 2 . In QCD the double- zero of the /3-function acts as a repulsive point for 
trajectories with a > and as an attractive one for a < 0. The domain in the complex a-plane 
covered by these trajectories will be called the "perturbative phase". Within this domain 1(a) is 
an analytic function but explicitly different from the physical answer (|2.9|). So, we define 



rPT/ \ f a da' a' [ f a da" ) ,„ . 

' "f^l^Kwf (2 ' 12) 



With the decrease of the phase of a, the radius (|2.11 ) increases and trajectories start to distort, 



being affected by the higher terms in the /3-function. At some critical point (r ~ 1, in general) a 
bifurcation occurs and the trajectories switch to another attractive point, different from —0. Such 
a jump results in a singularity (non-analyticity) of 1(a). It is easy to see from ( |2.11| ) that this 
happens for almost real initial a values: 



Ima c\a\ 2 



Re a Re a 



c\a\ w cRea < 1. (2.13) 



In the general case, the topological structure of the analyticity domains and, correspondingly, 
the number of jumps 1(a) experiences while approaching the real axis, can be quite complicated, 
reflecting the structure of the attractive zeroes of the /5-function in the complex a-plane. For the 
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Figure 1: Examples of RG trajectories (dashed) and separatrices (solid) for the polynomial ft- 
function (2.14a) with ai = oo, a = 1 (left) and for a singular ft (2.14b) with 7 = 0, h = 1 (right). 

sake of simplicity we shall restrict our wording to the scenario when the very first bifurcation ( "phase 
transition" from the perturbative domain) leads us directly to the physical phase; the generalization 
is straightforward. 

To illustrate this phenomenon one can consider two simplest examples with the two attractive 
points, —0 and a, with both examples respecting the first two PT terms of the QCD /^-function 
( p.4|) . In the first model with a polynomial ft the physical coupling freezes at a finite value a: 

ft(a) = — a 2 (l + a/ai)(l — a/a) , a^ 1 — 6T 1 = 7 . (2.14a) 

The second model employing a rational /3-function possesses a pair of complex conjugated poles 
and yields a(k 2 ) increasing logarithmically with k 2 — >0, that is a = 00: 

a 2 

P(<*) = -7 1 ^ tt-^ ■ (2.14b) 

Hy ! (1 - |7a) 2 + h 2 a 2 K ' 

In Fig.l examples of trajectories are shown together with the characteristic lines - separatrices 
- which border PT and PH domains. Crossing a separatrix in a causes non-analyticity in 1(a). It 
is worth reminding that the fact that the separatrices emerge from the origin +0 along the circular 
arcs, is of the most general nature and does not depend on the details of the model chosen for 
illustration. Thus, we arrive at the main conclusion of this Section, that 1(a) (for a rather trivial, 
purely perturbative, reason) cannot be analytic in any sector with a finite opening angle 8q covering 
the positive real direction in the a-plane. As a result, the usual Borel representation of the physical 
quantity J PH (a) does not exist. Indeed, if the integral 



1(a) = / duB(u)exp< > 

Jo I a J 



existed for some a = ao, it would define the function analytic for all a with Re a -1 > a^ 1 which is 
the interior of the circle |a — |ao| = |a . 

This non-analyticity in 1(a) at arbitrarily small positive a manifests itself in a singularity of the 
Borel image on the real positive w-axis, which is equivalent to the same-sign factorial asymptote 
of the PT coefficients. In the following we show that the Borel resummation actually yields the 
perturbative function (2.12), 1(a) = I PT (a). 
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2.2 Borel transform for the PT phase 

Following the above line of reasoning, one concludes that the analyticity domain in the perturbative 
phase, contrary to the physical phase, is broad enough to allow one to represent J PT (a) in terms of 
the Borel integral. Such a representation, however, involves imaginary values of the Borel parameter 
u. Namely, for sufficiently small a in the upper half-plane (lma>0) one can write 



J duB(u) exp|--| . (2.15a) 



The inverse relation gives the Borel transform B(u) as a regular function (for imaginary u), 

/oo— i/2r A<v 
^— I^{z~ l ) exp{uz} . (2.15b) 
-oo-j/2r 2m 

The integration line in z = 1/at shown in ( |2.15b| ) corresponds to a (small) circle in the upper 
half-plane of a, namely, \a — ir\ = r. 

In terms of the PT coefficients, one substitutes the (asymptotic) PT expansion for /, 



rPT _ 

n=l 



(2.16) 



into ( |2.15b| ). Closing the contour by the (multiple) pole at z = (Reuz < for Imz > 0), one 
obtains the standard improved series for the function B(u) 



n=0 



B( u ) = ^^±u n . (2.17) 



Thus, for imaginary values of u this series determines the regular Borel image B{u). 



2.3 "Condensate" contribution 

The IR renormalon problem arises when one attempts to reconstruct the observable 1(a) for real 
a by means of the perturbative series ( |2.16| ) processed via the Borel machinery ( |2.17| ) and the 



integral representation ( p. 15a ). The latter now runs, however, along the real positive w-axis. Then 
B[u) shows up a singularity (a pole, for the one- loop /^-function, or a cut in general) at u = Uq 
(in our normalization for the coupling, uo=p). As a result, the answer for / becomes ambiguous 
and, generally, complex, depending on the way one choses to pass by the singular point on the 
integration line. 

The PT coefficients a n do not care whether the expansion parameter a has been chosen real or 
imaginary. Therefore, we come to conclude that the IR renormalon is nothing but the outcome of 
an unjustified attempt to force the perturbative answer 7 PT ( |2.12| ) analytically continued outside 
its native phase, to represent the true 7 PH . Now we are in a position to cure the wound. To do 
so we observe that the physical answer as given by ( |2.9| ) can be reconstructed as a sum of two 
contributions: the PT piece, 7 PT (a), analytically continued to real a values plus the "condensate" 
(or "confinement") contribution originating from the discontinuity due to crossing the separatrix 
on the way to the real axis. Indeed, let us split the original integral along the physical trajectory 



in (|2.9| ), running from a to a into two pieces, — — ► a and a — > —0, to write 

J PH («) = J PT («) + H{a) , (2.18) 
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where 

^ da' a' f r a da" ) „ n , s r° da' a' 



The "condensate" contribution H satisfies the homogeneous RG equation ( |2.6|) . Therefore, turning 
from the a representation to the momentum dependence, one observes a pure power 

H(a) = CF{a) = |L, (2.20) 

with a complex constant C (or, equivalently, a dimensionful constant C\) to be determined from 
(|2.19 ). For example, for two models fl2.14| ) one has 



n^^-L}^™*^™ ( , 21a) 

and 

/(«) = a-7 exp |_I + (ll + h ^j a y (2.21b) 

The "condensate" constant C is expressed then, respectively, in terms of the hypergeometric or the 
Bessel function, depending on p and the parameters of the /3-function §1 . 

From ( |2.20| ) the reason becomes clear why do we refer to H as the condensate contribution 
(without quotation marks, from now on). Let us draw the reader's attention to the fact that J PH 
and J PT have identical PT expansions, and H has obviously none. 



3 Rescuing Borel representation: mission impossible 

A question arises, whether the Borel-type integral representation for the physical answer can be 
rescued. At least within our simplified approach in which the NP dynamics has been embodied 
into the behaviour of the /3-function at a ~ 1, the answer to this question is: in principle, yes. 
However, as a more extensive analysis shows, the Borel construction one is looking for proves to be 
not universal, depending essentially on the NP features of the theory. In other words, getting hold 
of the best possible PT information appears to be insufficient for this purpose, calling for genuinely 
new, non-perturbative, input. The problem is mathematically more involved, so that in the present 
note we give only a brief sketch of the procedure. 

One starts by explicitly constructing the Borel image of J PT . This can be done by translating 
the RG equation (|2.2|) into an equation for B(u). To this end one writes Q2.15b|) in the form 

B{u) = — / *L e **I PT (z- 1 )4-e {v - p)z (3-1) 
u — p Jc 2%i dz 



and, upon integrating by parts, makes use of (|2.2|) to arrive at the following compact symbolic 
equation 

uB(u)+(p(J^{(p-u)B(u)-l} = 0; d u = -^. (3.2a) 

The function here quantifies the deviation of the /3-function from its one-loop expression and is 
given by 

^ = r< YZ 2 ■ ( 3 - 2b ) 
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For the case of the pure one-loop /3(a) = —a 2 , one has 1/0 = 0, and one gets a simple pole solution 



1 1 / u \ _1 / Aixp \ 

B(u) = = — I 1 I , Uq = p \ = ——- in the standard normalization . (3.3) 

p-u p\ u J \ fa J 



Given a (rational) /^-function, it is straightforward to reduce (|3.2a ) to the differential equation (with 



proper initial conditions) for B(u). It is important to emphasize that the master (integral) equation 
(p.2|), as well as its differential counterpart, is an equation with the only singular point 



u 

x = 1 = 0, 

which, in general, results in a branch point in the solution, B(u), at u = u . Near the singularity 

B(u) = x^Fix) « c 1 , (3.4) 

V uqJ 



where 7 is given by the first two loops of the /3-function, see ( |2.4j ). Thus, one observes that both the 
position and the nature of the Borel singularity have purely perturbative origin. It is this universal 
singularity ( |3.4|) that governs the same-sign factorial growth of the PT series. 

For example, within the models ( |2.14|) one arrives at a second order differential equation. In the 



first case ( |2.14a| ) it is a homogeneous confluent hypergeometric equation for the function F(x) = xB, 
with the initial conditions F(l) = 1, F'(l) = 0. For the second model ( |2.14b|) B(u) is given by the 
solution of the inhomogeneous Bessel equation (see M for details). 
Having obtained B(u), one then has to construct the Borel integral 



rPT («) 



/ du B{u) exp I — 1 , (3.5) 
Jo la J 



passing the singular point, say, from above and to add the condensate contribution due to crossing 
the upper half-plane separatrix. From the Borel plane point of view, the latter contribution, H(a), 
as a solution of the homogeneous RG equation ( |2.6|) , can be written (up to a factor) as an integral 
of the same function B(u) along a quasi-closed contour, namely, around the cut uq < u < 00. Thus, 
the full answer can be represented as a sum of two integrals, one from zero to plus infinity (|3.5| ) and 
the second embracing the cut ( |3.4| ). The imaginary parts of these two contributions clearly cancel, 
as together they constitute the physical answer J PH . 

How much of the real stuff remains? To address this problem we restrict ourselves to a polyno- 
mial /3-function of power n+2. In this case one has the n-th order differential equation to determine 
B(u) (generalized confluent hypergeometric equation), with the initial conditions 



= 0, k= 1. ..n- 1. 

u=0 



pB(0) = l, (J^J (p-u)B(u 

Let us examine the large u asymptote of the solution. For u — > 00 equation ( |3.2a| ) reduces to 

1-0(O PidZ l )uB(u) = , 





which implies 



uB(u) ~ d Q exp j~r j + X! ^ ex P {~~~} ' fi-fy 



Here 7^ is a set of zeroes of the /^-function, other than the relevant one, the physical fixed 
point a. This particular term should cancel, however, in the combination of the ordinary (PT) and 
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the contour (NP) integrals (the discontinuity of B(u) is a solution of the homogeneous equation as 
well). Otherwise, the answer would be singular at a — a, which is not the case, since in the vicinity 
of the fixed point it is analytic: 

/ PH (a) = - + O(a-a) . 
V 

Moreover, for quite a while 1(a) stays regular above the fixed point, a > a, before the RG trajec- 
tories betray a for another, higher attractive point, i.e. another (unphysical) NP phase occurs in 
which the coupling is never small (no asymptotic freedom, that is). 

From this consideration one concludes that if a were the only zero, the condensate contribution 
(contour integral) would have to cancel the PT integrand above the singular point uq completely. 
This is simply because there would be no other terms in the sum ( |3.6| ) to maintain the asymptote. 

This explains a miracle of the finite Borel representation one obtains within the model for the 
two-loop /^-function with an (anti-QCD) negative 7. One can get this model from ( |2.14a| ) setting 
«i = 00 (7 = — a -1 ). The master equation (|3.2| ) for the Borel image then becomes the first order 
differential equation, 

(j)(d u ) = — d u = a d u ; a— [ (p - u)B(u) ] + uB(u) = , B(0) = 1/p . (3.7a) 
—7 du 



Its solution, 



B(u) = -(l--) exp - =- !-_ exp - , (3.7b) 

p \ pi la) p \ u J la: 



reveals explicitly the general features, namely, the nature of the singularity at u = p and the 
expected exponential behaviour at infinity. 

To evaluate the condensate contribution one first puts l/a\ = in fl2.21a|) to fix 



f a-a y 1 /* f 1 

f(a) = exp <^ 

a / la 



Then, one performs integration in ( |2.19| ) to obtain 
C = - 



* da' a' r a da' ( a - a'Y' p J p 
P{a') /P(a') ~ 7-o a'(l - a' /a) \ a' ) ^la' 



C dz (z - l)- 1+ > exp{ lp z} = (-Y P T( lp ) e— > ; 7p = -p 7 = £ > 

J-oo-iO \lpj a 



(3- 



Finally, 

H(a)=[-\ r( 7p )e->.(— — j -expj-^k (3.9) 

This expression, upon inspection, differs only by sign from the part of the ordinary (PT) Borel 
integral from Uq to 00, 



H(a) = - 



00+io 



duB{u) exp j j . 



As a result, for the physical answer in this model the finite-support Borel-type representation holds, 

J PH («) = £°duB(u) expj-^j. (3.10) 
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This observation has been independently made in ||. 

We now try to examine on this explicit example the standard routine employed in the renormalon 
analysis. Namely, it is accustomed to take the PV of the standard Borel integrals which amounts, 
for this model, to taking the value 



I B (a) = I PH (a) -ReH(a) = I pu (a) - f p (a) f-^-J r( 7p ) cos7r 7p . (3.11) 

At the same time, the uncertainty associated with the Borel resummation procedure is usually 
estimated as 

1 



51(a) = — 

ZTT 



duB(u)e- u/a 



c 



(3.12) 



with the contour running around the cut, which in our case yields 



/ e \ 7p 1 

61(a) = F(a) - Y( lp ) - |shi7r 7p | . (3.13) 
\lpj n 

Both the error of the PV Borel summation Ab = I pr —Ib and the estimated uncertainty 81 have the 
same correct power behavior (A 2 /Q 2 ) p . However, their relative magnitude, in general, mismatches: 

^- = r^—^ (3-14) 

ol |tan7T7 p | 

so that the actual error can become much larger than the estimated one, which happens, for example, 
when 7 P is numerically smallQ. The origin of such a mismatch is readily understood o. 

It is worthwhile to mention another interesting feature of an interplay between PT and NP 
contributions. Namely, when 7 p happens to be a positive integer, the singularity fl3.4|) disappears 
and B(u) remains analytic in the entire complex w-plane. There is no ambiguity in chosing the 
path of integration, and the Borel summation yields unambiguous - but incorrect - result! What 
is more intriguing from the viewpoint of the naive PT analysis is that the factorial growth of the 
coefficients in the expansion of 1(a) disappears: 

a n ~ n^- l (a)- n . (3.15) 

Resummation of the "subleading" 1/n corrections to the asymptote of a n kills factorials altogether. 
Analyzing mere perturbative series one would not infer any deficiency of the expansion, although 
the presence of the "NP condensate" has been demonstrated explicitly. 
For illustration, in the simplest case 7 P = 1 one has 

^-^^-M 1 -!})- (3 - 16) 

The first term is what one gets by the PT summation whereas the second one is the NP contribution. 
The former is merely a geometric series in a, showing no indication of NP effects. The latter term, 
on the contrary, has literally no PT expansion. In general, however, the situation is different and 
such an apparent splitting is absent, so that PT and NP pieces cannot be explicitly separated. 
It is just the genuine case of intrinsically mixed PT and NP contributions when the PT series 
normally send the message, via IR renormalons, about the presence of the non-perturbative effects. 

1 We parenthetically note that it is actually the case for the limit /3q — > oo (or nt — > — oo), viz., one gets 7 ~ l/rif 
if assumes that the higher order terms in the /?- function can be in turn obtained by only the leading in 
contributions, though literally one gets positive 7. 
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Notice, that for a generic /^-function (n> 1) an integer 7 does not mean analyticity: Mo becomes a 
logarithmic branch point of B{u). 

A curious observation is that for a given (polynomial) /3-function one can construct a polynomial 
P n (a{k 2 )) to replace a{k 2 ) in the definition of the "observable" Q2.1J ), such that B{u) becomes 



analytic and, thus, the PT expansion factorial-free H. However, such a construction proves to be 
non-universal with respect to p. 

In conclusion, let us mention that an attractive scenario of the finite-support Borel representation 
yielding I pu (a) analytic everywhere but at a = 0, is a peculiar property of the oversimplified two- 
loop model. Elsewhere such a possibility may accidentally occur only at specific, contrived values 
of the parameters. However, this would hold, once again, for one particular value of p, while for 
observables with the canonical dimension other than 2p the contribution from u>uq gets resurrected. 
Therefore we consider such an accident carrying no physical significance. 



4 Toy OPE for the Toy model 

For short distance Euclidean observables the Wilson OPE is known to provide a proper framework 
for describing power suppressed effects. It automatically cures, once and forever, the IR renormalon 
trouble for the price of modifying perturbative coefficients ||. The latter are changed only a little 
in low orders but become completely different at n—>-oo. How does it happen within our toy model? 

The analog of the OPE for 1(a) is introducing a normalization point fi 2 <C Q 2 (but still belonging 
to the small coupling domain) and considering 

rid/ \ T sd{ \ _ dk 2 f k 2 \ 2 rQ dk 2 ( k 2 \ . 2 \ 



1(a) = />) + />) - Jf w «M +1, w {q 2 ) ■ (4.D 

Now one treats perturbatively only the I sd (a) part. Following the above analysis, we want to study 
analytic properties of I sd (a) at a = 0. It is important that /i is kept fixed so that a(/i 2 ) depends on 
a via 

>(m 2 ) da' Q 2 

It implies that the RG trajectories employed for calculating I sd (a) are always of the fixed finite 
time interval r = logQ 2 //i 2 . This is what makes the crucial difference: I sd (a) is an analytic function 
at a = with a finite radius of convergence! 

In reality, one is interested in a proper value of /i necessary for the OPE which, for practical 
reasons, is desired to be as low as possible. One should bear in mind that now the coefficients of the 
perturbative series for I^f depend explicitly on the ratio fi 2 /Q 2 . Under such circumstances a critical 
momentum scale /z sri emerges [0] such that the series converge for a<a(ii 2 ) if \i is chosen above /i sri , 
and diverge (for sufficiently small a, that is, large \n(Q 2 / fi 2 ), sid) if /i</i s h- We denote the value of 
the coupling corresponding to this scale by a sri = a(/i s h) < «• Leaving aside tiny details concerning 
perverted /5-functions H], one can determine this OPE-b order line value solving the equation 



sh -p(a') 



where the integration contour travels to infinity along a separatrix. Such an equation always has 
a unique solution between zero and a. In a general case when the structure of zeroes is rich and 
there is a set of separatrices, one should take the minimal solution. (Singularities of the /^-function, 
if any, should also be considered, with the position of the singularity replacing 00 in the integral in 
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one obtains a s ^ ~ 0.782 a. For the model ( 2.14b ) 



( |4.3|) ). For example, in the two- loop model 
with 7 = 0, discussed in 0, a s h = h~ x . 

Borel non-summability in I within the OPE lies in the second, large distance, piece I ' , which, 
of course, cannot be even dreamed to be calculated in terms of a expansion. This contribution is 
explicitly associated with the particular, physical, phase. 

We note that an attempt to do the OPE without an IR cutoff, i.e. putting fi—>-0 (or, equivalently, 
using Dimensional Regularization to "renormalize" the power convergent integrals like /) means 
calculating I PT (a) instead of I^(a), and leaving only H(a) in i^ d (a). In other words, it is an attempt 
to entirely subtract "perturbative corrections" from 7 ld . It is just this action that generates the IR 
renormalons in the coefficient functions making them uncalculable. As long as the exact analytic 
expression for 1(a) is known, such a routine seemingly does not pose particular problems. However, 
in practical applications it looks rather dangerous when multiloop corrections are incorporated. 
Indeed, consider for example the Q 2 -dependence of the observable. Translating the behaviour in 
a(Q 2 ) into the Q 2 -dependence one finds a smooth behaviour of I at small Q 2 : 



iHQ 2 )) 



a: 



P 



1 



dj3(a) 



P + P 



A 2 

, QCD, 



da 



> ; for Q < A QCD . (4.4) 



On the contrary, according to ( |2.20| ), both the "perturbative" and "purely non-perturbative" parts, 
taken separately, are singular: 



J PT (a(Q 2 )) ~ H(a(Q 2 )) oc 



(4.5) 



Therefore, the naive Borel resummation typically yields a strongly corrupted picture of the actual 
scale of long distance phenomena. We will return to this important point in [[|. 



5 Conclusions 

We argued that the IR renormalons have a transparent mathematical origin and are associated with 
employing the same short distance coupling as an expansion parameter for describing physics at 
both large and small scales. They do not actually show what happens at low momentum scales 
but only send a message that something may happen there. Inspired by this idea » we studied a 
simplified model for hard QCD observables in which all the details of strong dynamics are embodied 
into an IR non-trivial /^-function. 

We explicitly showed that, contrary to a rather popular believe, IR renormalons are not directly 
related to the Landau singularity in the running coupling. Even if it freezes at a small value, the 
same-sign factorial growth of the PT coefficients remains intact. At the same time, there are curious 
examples when no IR renormalon is around (the PT series converges!) but the PT expansion is as 
deficient as in a general case. 

The position and the nature of the IR renormalon singularities are determined by the first two 
coefficients of the /3-function, i.e. are governed by the deep PT domain. The true origin of the IR 
renormalons is rooted in a non-trivial phase structure of QCD. 

We examined how well can renormalons represent, in general, the actual long distance effects 
and found them unsatisfactory in many important respects including the absolute magnitude, an 
estimate of incompleteness of the purely PT approximation and the steepness of the Q 2 -dependence 
in the transition regime from short distances to the strong interaction domain. 
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As far as the ultraviolet renormalons are concerned, in the framework of our model they appear 
in the observables represented by an integral 



Straightforward application of the analysis outlined in the present paper shows that U is analytic 
at small a except for only a narrow beak around the real negative axis. Therefore, its Borel 
resummation yields the correct result. This being a positive statement, we are not sure of how 
heavily does the conclusion rely on particular features of the oversimplified model. Therefore we 
refrain from definite claims of its applicability to actual QCD. 

When this paper was in writing we learned about the preprint by G. Grunberg || where the 
finite support Borel representation has been constructed for the two-loop ("non-QCD") /3-function, 
in agreement with our finding. 
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